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ABSTRACT 



Aims. The so-called Limber equation is widely used in the literature to relate the projected angular clustering of galaxies to the 
spatial clustering of galaxies in an approximate way. This note gives estimates of where the regime of applicability of Limber's 
equation stops. 

Methods. This paper revisits Limber's equation, summarises its underlying assumptions, and compares its predictions to the 
accurate integral of which the Limber equation is an approximation for some realistic situations. Thereby, the accuracy of the 
approximation is quantified. Diff'erent cases of spatial correlation functions are considered. 

Results. Limber's equation is accurate for small galaxy separations but breaks down beyond a certain separation that depends 
mainly on the ratio a/vm and to some degree on the power- law index of spatial clustering; a is the la- width of the galaxy 
distribution in comoving distance, and Vm the mean comoving distance. Limber's equation introduces a 10% systematic error 
that may be important at a scale of about a few degrees or in some cases, narrow galaxy distributions for example, even at 
arcmin scales. As rule-of-thumb, a 10% relative error is reached at ~ 260 a /rm arcmin for 7 ^ 1.6, if the spatial clustering is a 
power- law ^ cx r More realistic ^ are discussed in the paper. 

Conclusions. Limber's equation becomes increasingly inaccurate for larger angular separations. Ignoring this effect and blindly 
applying Limber's equation can possibly bias results for the inferred spatial correlation. It is suggested to use in cases of doubt, 
or maybe even in general, the exact equation that can easily be integrated numerically in the form given in the paper. This 
would actually eliminate all problems discussed in this paper. 

Key words, galaxies - statistics : cosmology - large-scale structure of Universe : cosmology - theory 



1. Introduction 

If the distance of individual galaxies is not known, but only 
statistically for a sample of galaxies on the whole, one has 
to infer the real-space correlation from the observed angu- 
lar correlation on the sky. This task has been performed 
many times since the early days of g alaxy surveys as, 
for ex a mple, in the pioneering work of Totsuii fc Kihara 
(196i), iGroth fc Feeble^ (|l977l) and iDavis et all (|l977^ 
based on the Lick galaxy survey. To relate the pro- 
jected correlation function, lu{6), to the spatial, three- 
dimensional correlation function, ^, man y author s use an 
approximation that was introduced by iLimbej ([l953). 
This approximation is nowadays referred to as Limber's 
equation. 

It will be shown in Sect. I^T^ that Limber's approxima- 
tion diverges when galaxies are closely distributed. That 
Limber's equation breaks down for narrow distributions is 



* Founded by merging of the Institut fur Astrophysik 
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no surprise, as it is assumed that the distributions do no 
vary considerably over the coherence length of the spatial 
correlation function. But how wide has a distribution to 
be? How much do we bias our estimate for the angular 
correlation by using Limber's equation, and is there an 
alternative? These questions are addressed by this paper. 

Peebles! (Il980l p.l99) argued that this break-down oc- 
curs at separations of roughly 30°, but made assumptions 
about galaxy surveys that do not meet the standards of 
many contemporary wide-field surveys. Due to the avail- 
ability of spectroscopic or photometric rcdshifts quite nar- 
rowly binned distributions arc often used in current stud- 
ies. 

This paper delineates the circumstances in which 
Limber's approximation is valid and gives an alternative 
for situations where it is not valid. Throughout this paper 
a flat universe will be assumed. 

The exact equation for lo is given in Sect. 12.11 The 
thereafter following Sect . 12 . 2l summarises Limber's approx- 
imation and investigates why it breaks down for peaked 
galaxy distributions or beyond a certain angular separa- 
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tion. Alternative approximations for cases of narrow dis- 
tributions of galaxies or for larger 9 are given inside the 
Sect. 12.31 By looking at different distributions with vary- 
ing widths and means, the accuracy of the Limber ap- 
proximation is assessed in Sect.|3l The paper finishes with 
conclusions in Sect.Hl 

2. Relation between spatial and angular 
correlation function 

2.1. Exact relation with only some restrictions 

Consider two number density fields of galaxies, ni(r) and 
712(7"), as a function of the position r within the comov- 
ing frame. One defines the spatial clustering, or spatial 
correlation function, ^(Ax) in terms of the average num- 
ber of pairs of ni- and n2-galaxies that can be found 
when considering small volumes dVij^ with a separation 
of Ax = \ri - ral: 

(ni(ri)n2(r2))dVidV^2 = (m) (n2> [1 + e(Ax)]dl^idF2 • (1) 

(.) denotes the ensemble average. It is already assumed 
here that the random fields ni.2 are statistically homo- 
geneous and isotropic so that ^ is a function of the dis- 
tance Ax only. If we have ni — n2, then ^(Ax) is a auto- 
correlation function of a galaxy number density field. 

Observing the galaxy density fields as projections on 
the sky, 



MO) 



drqi{r) {r6,r) 



(2) 



yields a different correlation function, w{6), when pairs of 
galaxies, fii and n2, within solid angles dfli/2 of angular 
separation 9 are considered: 

(ni(6'i)n2(6'2))dr!idr!2 = {fii) {n2) [I + Lj{e)]dnidn2 ■ (3) 

Here, qi{r) is a filter or weight function (of the number 
density) solely depending on the comoving radial distance 
r. In the context of galaxy number densities, this filter 
needs to be normalised to one, i.e. 



drqi{r) = 1 



(4) 



it specifies the frequency of galaxies within the distance 
interval [r, r -I- dr[. The details of the filter depend on the 
selection function of the survey. The vector 9 is the di- 
rection of the line-of-sight. The observer is sitting at the 
origin of the reference frame, r = 0, while other points 
on the line-of-sight are parametrised by r9. To take into 
account a possible time-evolution of rii the second argu- 
ment in ni{r6, r) is supposed to parametrise the look-back 
time, t{r), as a function of the comoving radial distance 
at which the random field is observed: 



t{r) ~ - dr' a(r') , 
c Jo 



(5) 



where a(r) is the scale factor at comoving distance r, and 
c the vacuum speed of light. 



As already mentioned the arguments stated above as- 
sume a Euclidean geometry for the Universe. This is not 
a strong restriction, though, because any other fiducial 
cosniological model can be incorporated by redefining the 
filter Qi - b y absorbing F (x) in the (relativistic) Limber 
equation of IPeebleal (llQSnt l (Chapter III, Sect. 56) into qt 
- which leaves the form of ([2]) effectively unchanged. 

Let us now derive a relation between lo and ^. Imagine 
we observe the spatial fluctuations ({ni){r) is the mean 
density at cosmic time t(r)), 



0i,2\r6, r) = — 1 



(6) 



(ni,2)(r) 

of the two random fields as projection on the sky towards 
the direction 6: 

/>oo 

hAO) = / drpi,2(r)5i,2(r0,r) . (7) 
Jo 

Here, pi are the filters of the spatial fluctuations, ^i, that 
correspond to qi (the filter of the density fields) \2 



Pt{r) = qt{r) {n^){r) 



drqi{r) {ni){r) 



(8) 



This relation between qi and pi follows from the definition 
of the density contrast of the projected (number) density: 



1 



(9) 



If the mean density is constant, both associated filters pi 
and qi will be identical. 

According to the definition of to we have for sight- lines 
i?i and i?2 spanning an angle 9 — <i9i, '02- 

u{e) = (5i(i9i)(52(i52)) (10) 



dri I dr2Pi{ri)p2{r2){5i{ri-di,ri)52{r2-d2,r2)) 
Jo 

dri / dr2Pi{ri)p2{r2)^[R, j , (11) 



^0 



where 



2rir2 cos6' 



(12) 



As before the second argument of the spatial correlation, 
^{R,r), parametrises the time, t{r), at which the spatial 
correlation function is observed. 

Two assumptions had to be made to arrive at (|lip : a) 
the random fields (5i^2, and hence also their projections, 
Si^2, are statistically isotropic and homogeneous, and b) 
the time-evolution of ^ is small within the region where the 
product Pi{ri)p2{r2) is non- vanishing. Due to assumption 
a) Lo depends only on 9 and is independent of the directions 
1^1,2- Owing to assumption b) we can approximate the 
spatial correlation of fluctuations at different cosmic times 
(radial distances) t{ri) and t{r2) by a representative ^ at 
time (ri -|- r2)/2. 

Note that Eq. (fTTj) is, under the previously stated as- 
sumptions, valid even for large 9. 

^ In this paper the focus will be on the angular clustering 
and hence normalised pi are assumed, albeit all equations and 
the conclusions are also valid for other types of projections. 
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2.2. Limber's approximation and its breakdown 
Limbej (|l953l) introduced several approximations in order 



to simplify the integral Eq. pT|) . The details of this ap- 
proximation are given in Appendix[Xl Here only the result 
is quoted: 



ui{9) = / drpi{r)p2{r) 



+ 00 



dAr ^{R,f) 



R = \/f'^0^ + Ar2 



(13) 
(14) 



The useful Eq. has frequently been used in the 
astronomy community because it al lows one to find an 
analytical expression ()Peeblesl \l98dt) for uj in case of a 
power-law like 



[vq is the clustering length) namely 

1-7 



IRAD 



The angular clustering amplitude at 6* = 1 RAD is 



(15) 



(16) 



- r(7/2-l/2) 

r(7/2) 



dr pi{r)p2{r)r^ . (17) 



It has been found that a power-law ^ is a fairly 
good description of the tr ue spatial correlation function 



(e.g. IZehavi et al. 20021 ) of galaxies for scales below 



15Mpc/h, which makes Limber's equation quite prac- 
tical. Typical values are 7 ^ 1.8 and rp ^ 5/i^^Mpc. 

However, due to the integral, Ai^ diverges to infinity if 
the width of the distributions pi,2 approaches zero and as 
long as pi,2 remain overlapping; the latter is trivially true 
if we are considering auto-correlations, pi = P2- To make 
this point, let us assume that pi^2 are top-hat functions 
with centre rc and width 2Ar: 



1 



2Ar 



1 for r G [rc 
otherwise 



Ar, Tc + Ar] 



It follows for A,, 



4A7" 



rc+Ar 



— Ar 



2Ar 



(18) 



(19) 



The last step is valid for small Ar <C rc, which shows that 
diverges for narrow distributions as Ar^^. This small 
calculation implies that Limber's equation possibly over- 
estimates the angular correlation oj to some degree. 

In the astronomical literature, Eq. (|13p is also known 
in other fo rms involving the three-dimensiona l power spec- 
trum (e.g. Kaiser 1998t Hamana et al.|[2004 ). No matter 



the form of Limber's equation, it always suffers from the 
divergence previously discussed and from the inaccuracies 
that are going to be discussed in the following. 



2.3. Thin-layer approximation 

Approximation Eq. (jlSp fails if the weight functions be- 
come delta-like functions, i.e. for galaxies being located 
inside one layer at comoving distance rc. To obtain the 
correct solution for 



Pi,2(?') = Snir - ri,2,c) 



(20) 



(where Suir) is the Dirac delta function) we simply have 
to go back to Eq. (fTT|) and find: 



'2,c 



2ri cr-2,c cos( 



r2x 



(21) 



For ri c = r2,c = ^'c, i-e. for auto-correlations which are 
mainly considered in this paper, we obtain as quite simple 
and intuitive result 



cos 6*, r^ 



e(^c 



(22) 



This means that for very narrow filter we essentially ob- 
serve a rescaled ^. An immediate consequence is that lj 
has the same slope, 7, as ^, if ^ happens to be a power-law 
as in Eq. (jl5|) : on the other hand for wide filter, in the 
regime where Limber's approximation is accurate, lo has 
a shallower slope of 7 — 1 . 

It can be shown, as done in Appendix IB| that the thin- 
layer approximation is not just a sole particularity of ex- 
tremely peaked weight functions p. The solution is asymp- 
totically approached by (fTTj) if the separation 6 gets large 
enough, namely 9 > \/2;r- hi the case of auto-correlations 
p{r) = pi{r) — P2{f')] o is the la-width and r^ the mean 
of the weight p{r) (both in units of comoving distance). 

This discussion shows that Limber's equation and the 
thin-layer equation are approximations for two extremes, 
namely for small 6 and large respectively. The transition 
point between the two regimes is mainly determined by 
the ratio cr/rm implying that the accuracy of Limber's 
equation is, in the case of auto-correlations at least, mainly 
determined by this ratio. 

2.4. Numerical integration of the exact equation 

With the availability of computers there is no longer a 
need to use any approximations in this context. However, 
one finds that solving Eq. pT|) for lA.2|l by numerical inte- 
gration with a power-law ^ is cumbersome because many 
fine bins are needed to achieve the desired accuracy. This is 
mainly due to the argument i? in ^ which does not sample 
f on a log-scale, which however should be aimed at with 
a power-law like 1^. It is advisable to consider a numerical 
integration with different coordinates such that the spatial 
part of ^ depends just on one integration variable R. 

This can be done by changing the integration variable 
in the inner integral of (IA.2|) from Ar to R: 



iu{9) = 



1 



cos ( 



df / dRQ{f,A)^{R,f) 



R 
A ' 



(23) 



f^2(l-cos9) 
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where 



and 



' 2f^{l - cos 9) 
1 + cos 6* 



(24) 



Q{f,A)=pi{f~A)p2{f + A)+pi{f + A)p2{f-A) .(25) 



The expression may look more complicated than 
(fTTj) or (|A.2p but is superior for numerical integrations: 
The integration over R can be done on a logarithmic scale, 
the integration limits for R can be adjusted for ^'s that 
effectively vanish somewhere in the integration interval. 
Furthermore, also the integration limits for f can be ad- 
justed if Q{f, A) vanishes effectively somewhere within the 
integration limits. 



3. Accuracy of Limber's equation 

The aim of this section is to assess the accuracy of 
Limber's equation by comparing its predictions to the 
exact solution for the angular correlation function. For 
this purpose, different spatial clustering models f are used 
that are constant in time. Moreover, the focus is on auto- 
correlations w, hence pi{r) — P2{r) = p{r). 

In order to quantify the accuracy of Limber's approxi- 
mation, the exact solution Eq. pTjl is solved by numerical 
integration, by means of (j23p , and compared to the result 
from Limber's equation (|13p . In order to make a fair com- 
parison, the same small angle approximation as in (fT3|) is 
applied to the exact equation, which as it is written in 
(j23|) does, of course, not require small galaxy separations: 



1 



1 -h cos 6* 



2 



(26) 

(27) 
(28) 



3.1. Model galaxy probability distribution 



The intention is to study the accuracy of Limber's equa- 
tion for relatively narrow distributions of galaxies. 



p{r) cx exp 



1 (r 
2 



(29) 



Thus, galaxies selected by this filter are centred on with 
some scatter cTp. The prefactor accounts for the appar- 
ent increase in the angular galaxy number density with 
distance r within the light-cone which only for moder- 
ately wide or wide distributions becomes relevant. The 
prefactor has been added to make the distribution more 
realistic. 



The parameter dp in (|29|) is not identical with the Icr- 
scatter of the distribution, and is not identical with Tc ; 
this is only roughly correct for small x = ap/vc. However, 




1 1 00 

angular separation in [arcmin] 



Fig. 1. Shown here is the angular correlation function 
Lj for three different ratios of the galaxy distribution, 
see Eq. (l29l) . Tc = 2/i~^Gpc. The spatial correlation, 
is a power-law with slope 7 = 1.77 and clustering ampli- 
tude To = 5.4/i~^Mpc as roughly found for real galaxies. 
The thick lines are the numerical integrations of the ex- 
act equation, Eq. while the thin solid lines are the 
predictions by Limber's equation, Eq. (fTB|) . The small an- 
gle approximations assumed in Limber's equation are also 
applied to the exact equation. The relative difference be- 
tween exact and Limber solution is plotted in the bottom 
small panel. The thin dashed line is the exact solution for 
an extremely peaked galaxy distribution, Eq. (j22p . 

it can be shown that the ratio y = a/r,^ is entirely a func- 
tion of X, namely approximately 



y 



O.Ola:^ -exp (-1.29 x^i-^^)] , (30) 



with an accuracy of a few percent for < a; < 2.5. This 
relation is used in the following to find a p{r) for given y 
and Tp. 



3.2. Power-law spatial correlation functions 

The angular clustering, lu, for three different galaxy dis- 
tributions has been computed for the Figure HI ^ is the 
same in all three cases. As expected. Limber's equation 
describes the exact solution quite well for small angular 
separations 6 but starts to deviate from the 7 — 1 power- 
law for larger separations, becoming steeper and tending 
to look very much like the thin-layer solution ([22l) (thin 
dashed line) which has the slope 7. Therefore, the exact 
angular clustering is a broken power-law if the spatial clus- 
tering is a power-law. This behaviour is expected from the 
discussion in Appendix IB] 

The break position is important for the Limber approx- 
imation: The closer one gets to the break position (and 
the farther beyond it), the more inaccurate is Limber's 
description. When ^ can be assumed to be a power-law 
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we can define a sensible position ^brcak for the break by 
working out where the Limber and the thin-layer approx- 
imation intersect. As shown in Appendix [Cl this is at: 

2^/3 



^brcak 



r(7/2) 



(31) 



r(7/2 - 1/2) rm 

This relation reflects what has already been seen in 
the foregoing discussion. The break position is mainly de- 
termined by the ratio a/r^ - but is also influenced by 
the slope 7 of ^: the prefactor of a/r^ in ([3T|) varies be- 
tween 0.31 and 1.18 for 7 — 1.2, 2.1, respectively: steeper ^ 
are projected more accurately by Limber's equation than 
shallower ^. 

The relative error of Limber's approximation at the 
break position is already large, roughly 100%. Therefore, 
the validity regime of this approximation stops well be- 
low the break position. As a rule-of-thumb, the relative 
error is about 10% at 6'broak/lO which has been estimated 
from several plots like Figure [T] Therefore, a 10% error is 
reached at about ~ 260 a /r^ arcmin for 7 ^ 1.6. 

3.3. More realistic galaxy clustering 

Although (spatial) galaxy clustering is empirically well de- 
scribed by a power-law on scales below r ~ 15 Mpc/h, ^(r) 
on large cosmological scales eventually adopts the shape of 
the dark mat ter clustering which d oes not have a constant 



slope 7 (e.g. Weinberg et al. 2004( l. Therefore, the previ 



ous discussion is only applicable if the effective physical 
scales, that correspond to the angular separation interval 
under investigation, are below r ^ 15 Mpc/h. In the oppo- 
site case, we expect that the accuracy of Limber's equation 
is not only a function of cr/rm and an average of 7, but 
also a function of rm- 

To obtain more realistic ^(r), two representative mod- 
els for galaxy clustering based on a halo-model as in 
Magliocchetti fc Porcianil (NFW parameters; 2003f) . are 
used in this section: one ^ is supposed to represent the 
clustering of galaxies that preferentially inhabit dark mat- 
ter haloes of larger mass (Mhaio > 10^^-^ Mgun] "red galax- 
ies"), another ^ describes galaxies that also live in smaller 
mass haloes (Mhaio > 10^^ Mgun; "blue galaxies"). These 
two ^'s are placed at different redshifts between < z < 5, 
are "observed" on the sky through filters p{r) of varying 
widths cr/rni, and compared to Limber's equation in or- 
der to work out the relative error. As to the accuracy of 
Limber's equation the amplitude of the correlation func- 
tions is irrelevant, it is only the shape that matters. For 
all redshifts considered the same ^ is assumed, that is the 
halo-model prediction for z = 0. To give an estimate of 
the relative error of narrow filters, this should be a fair 
assumption. See Figure [5] for a particular example of the 
angular clustering. 

Figure [3l shows the estimated separation limits. The 
dependence of the filter on the centre, r,,!, however can 
clearly be seen, becoming weaker for larger r^a- The differ- 
ence between red and blue galaxies is most pronounced for 
small r,„ but there are clear differences for a/rm < 0.02, 




angular separation in [arcmin] 

Fig. 2. Angular clustering (arbitrary scaling) of the red 
galaxy sample (upper panel, thick solid line) centred on 
z — 0.2, the relative width of the distribution in comoving 
distance is cr/r^ — 0.02. For comparison, also in the up- 
per panel, the clustering of the blue population is plotted 
as well (thin solid line). The dashed and dotted line are 
the approximate solutions according to Limber's equation 
and the thin-layer approximation (upper panel) . The lower 
panel depicts the relative error of the Limber (dashed) and 
thin-layer solution (dotted) for the red galaxies, respec- 
tively. 

even for large Tm- Limber's prediction of the angular cor- 
relations of red galaxies is more accurate than the predic- 
tions for blue galaxies because on megaparsec-scales ^ of 
red galaxies is steeper than ^ of blue galaxies. 

4. Summary and conclusions 

The Limber equation is an approximate solution for the 
angular correlation function of objects, mainly applied in 
the context of galaxies, for a given spatial correlation func- 
tion. 

In this paper, it was shown that the assumptions made 
for Limber's approximation become inaccurate beyond a 
certain angular separation, which essentially depends on 
(j/rm and the shape of the spatial correlation function. It 
should be emphasised that this inaccuracy has no relation 
to the flat-sky approximation (small separations) that is 
commonly adopted when working with Limber's equation. 

Beyond what separation 9 does Limber's approxima- 
tion become wrong by 10%? For power-law ^'s this sep- 
aration depends entirely on 7 and tr/rm of the filter (se- 
lection function) used. Spatial correlations with larger 7 
are treated more accurately by Limber's equation than ^ 
with shallower 7. As a rule-of-thumb, a systematic error 
of 10% is reached at about 260 — arcmin for 7 ~ 1.6. 

For more realistic ^'s with a running slope 7, the sep- 
aration is also dependent on r^, the centre of the filter. 
Even for relatively wide filters with a/r^n ~ 0.2, Limber's 
equation becomes inaccurate on the 10% level at a scale of 
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red galaxies 



0.01 0.02 0.05 0.1 0.2 0. 
relative width ot filter 



01 0.02 0.05 0.1 0.2 

relative width of filter 



Fig. 3. Galaxy separation, as a function of cr/rm of the 
filter, at which Limber's equation becomes inaccurate by 
10% for red (left panel) and blue (right panel) galaxies. 
Different lines correspond to different the correspond- 
ing redshifts are the line labels. Note that for estimating 
the relative error the local correlation function of red and 
blue galaxies was simply moved to higher redshifts. The 
(flat) fiducial cosmology assumes f^m = 0.3. 




Fig. 4. Filter (galaxy distributions) p{r) appear "wider" 
in redshift, Pz{z) = p{r{z))^, than in comoving distance, 
i.e. usually {O^ = cr^/z) > [9^ = c^/r^) where and 
are the la-widths in redshift, and comoving distance, re- 
spectively; is the mean of the weight in comoving dis- 
tance, while z is the mean in redshift. To give an estimate 
of the correction ^r/^z (y-axis) as function of 9^ (x-axis) 
Pz's are assumed that are close to Gaussians. Each solid 
line corresponds to different means z S [0.05, 7] (upper to 
lower line). 



a few degrees separation. On the other hand, for narrow 
filters of about a/r^ < 0.02, a systematic error of 10% is 
already reached after ~ 10' (or even a few arcmin depend- 
ing on r,n). 



galaxy sample 


2 ± (Tz 




Quadri et al. (2007) 






ii'-dropouts 


2 = 2 ... 3.5 


~ 0.16 


Lee et al. ('2006) 






fZ-dropouts 


3.2 ± 0.3 


~ 0.04 


B435-dropouts 


3.8 ±0.3 


~ 0.03 


Veoe-dropouts 


4.0 ± 0.3 


~ 0.02 


Overzier et al. ('2006') 






i775-dropout 


6.0 ± 0.5 


~ 0.03 


Adelbereer et al. (2005) 






BM 


1.69 ±0.30 


~ 0.12 


BX 


2.24 ±0.37 


~ 0.09 


LBGs 


2.94 ±0.30 


~ 0.05 


Hildebrandt et al. f2005) 






LBGs 


2.96 ±0.24 


~ 0.04 


Ouchi et al. ('2004) 






BRi - LBGs 


4.0 ±0.5 


~ 0.05 


Viz - LBGs 


4.7 ±0.5 


~ 0.04 


Riz - LBGs 


4.9 ±0.3 


~ 0.02 


LAE 


4.86 ±0.03 


~ 0.003 



Table 1 . List of typical values for the relative width of the 
survey selection function, a/r^ (comoving distance), in 
recent works using the dropout technique (LBGs: Lyman- 
break galaxies, LAE: Lya-emitter) . Note that typically 
only separations less than about one arcmin are considered 
so that the accuracy of Limber's equation is, apart from 
the only LAE study in the list, no real issue. 



Note that a filter can appear quite wide, i.e. a^/z is 
large, as a function of redshift z while it is actually quite 
narrow as a function of comoving distance, the latter be- 
ing the relevant parameter for the accuracy of Limber's 
equation. In particular, this means that applying Limber's 
equation to high redshift galaxies may be a problem. The 
reader might find Figure |4] useful to quickly look up the 
approximate conversion between a^/z and a^/r^. 

Table [1] shows a small list of some recent studies of 
the high-redshift universe using the dropout-technique. 
According to the table the relative width of the selec- 
tion filter frequently reaches values near 0.03. Therefore, 
Limber's equation should in these and similar studies not 
be used when separations larger than ~ 20' are considered. 

In cases where Limber's approximation is in doubt, 
that one uses the exact Eq. that can easily be inte- 
grated numerically. 

Limber's approximation is also used in gravitational 
lensing to quantify correlations between the gravity- 
induced distortions of faint galaxy ima ges (see e.g. 
iBartelmann fc Schneideill200ll: lKaiserlll998l ). In this case, 
the filter function relevant for lensing is the lensing effi- 
ciency that, quite naturally, has a wide distribution; one 
can roughly estimate cr/r,„ w 0.22 if the source galaxies 
have a typical value oi Zg — I. The centre of this filter is 
at about z w 0.5. These values will depend on the distri- 
bution of source galaxies in redshift and the fiducial cos- 
mological model. Taking these values as a rough estimate 
of the relative lensing filter width and centre, and assum- 
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ing that the dark matter clustering is somewhere between 
the clustering of red and blue galaxies, one can infer from 
Figure [3] that a two-point auto- correlation function of the 
convergence, based on Limber's equation, is accurate to 
about 10% for separations less than several degrees; be- 
yond that, an alternative description should be used. 
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Appendix A: Details of Limber's approximation 

In order to find an approximation of Eq. (jlip . one intro- 
duces for convenience new coordinates. 



ri + 7-2 



Ar 



r2 - ri 



(A.l) 



This change of coordinates renders (fTTj) in the following 
manner: 



o 

df 



+2r 



'0 J -2f 

where 



dArpi [r — 



P2 



^)e(i?',f)(A.2) 



R' 



Ar2 

2f2(l - cos 6*) + —(1 



(A.3) 



This equation is simplified further by making the following 
approximations which are roughly satisfied for wide weight 
functions pi.2, and for ^ that fall off sufficiently fast over 
the typical width of pi_2- 



/+2r r+oo 
dAr^{R',f) w / dAr^{R' 
-2f J-oc 



(A.5) 



Especially the approximation (|A.4p is characteristic for 
Limber's equation. It is justifie d if the weight functions 
Pi,2 d o not "vary appreciably" ([Bartelmann fc Schneider 



2001 



p. 43) over the coherence length of structures de- 
scribed by ^ ~ typically a few hundred Mpc in the context 
of cosmological large-scale structure which means we 
consider cases in which the coherence length is small com- 
pared to the width of the weight functions pi,2- In total, 
these assumptions lead to the (r elativistic) Limber equa- 
tion (|Limbeilll95i |Peebleslll9"80l ): 



where 



/>00 /' + 00 

f dfpi{f)p2{f) I dAr^{R,f), 

J -oo 



R=Vf' 



(A.6) 



(A.7) 



For historical reasons, as a further approximation it is as- 
sumed in the above equations that we are dealing with 
small angles of separation, 9, by which we can introduce: 



1 -I- cos 6* : 



1 — cos 9 



(A.8) 



which are the mean radial comoving distance and differ- 
ence of radial distances, respectively, of a pair of galaxies, we can work out what happens for larger 9. Note that 



These two approximations are accurate to about 10% 
for angles smaller than 9 < 40° which covers the typical 
range of investigated separations. Usually, when employ- 
ing Limber's equation this approximation is automatically 
used. 



Appendix B: Approximation for larger separations 

It will be shown here that the thin-layer solution of the 
Sect. 12.3] is in fact asymptotically approached if the sepa- 
ration 9 only gets large enough. 

To see why the exact solution behaves like this we have 
to go back to Eq. (jA.2p and in particular (|A.3p . From that 
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the following focuses on auto-correlations. The Ar-ternr 
in (|A.3p can be neglected compared to the f-term if 



A?" 
2^ 



■ tan 



(B.l) 
(B.2) 



Let us say the distribution p of galaxies has a char- 
acteristic width (variance) of a and a mean of r^. 
Pairs of galaxies from this distribution will typically 
have r = {{ri + r2)/2) K and for their mutual dis- 
tance Ar « {{ri — r2y) fa ^/2a. Therefore, the condi- 
tion (|B.2p will be given for 9 large enough to fulfil 



(B.3) 



or for small 9 (10% accuracy for 9 < 60°) approximately 
9 > V2— . (B.4) 



Thus, this is a reasonable approximation to use. 

The relation (jC.2p is useful because it allows one to 
estimate the expected position of the break in the power- 
law behaviour of lo. As the accuracy of Limber's equation 
is primarily an issue for narrow p we can derive from (jC.2|) 
a rule-of-thumb for the break position. For narrow p, it is 
sensible to assume a top-hat function as in Eq. (|18p . Then, 
employing (jC.SP and (jC.2p one gets: 



^brcak 



2^/3 r(7/2) o_ 
^ r(7/2- 1/2) ^ 



(C.4) 



Note that the variance of the top-hat 
a = Ar/-v/3, where Ar is the width of the top- hat. 



IS 



In this regime, where the Ar-term in B! is negligible, Eq. 
(jA.2p simplifies to 



uji9)=J dfF{f)^[rV2Vl-cos9, 

which is a family of thin-layer solutions 
eraged with the kernel 

/■+^' /_ Ar\ /_ Ar 

F{r) ^ J dArp\r + — j plr - — 



(B.5) 
weighted av- 

(B.6) 



Already for moderately wide p{r), F{r) is relatively 
peaked, so that for large enough 9, Eq. (jB.sp is close to the 
solution ([221) with Tc « r,n. Note that F{r) is normalised 
to one. 



Appendix C: Break position 

Limber's equation is an accurate description for small 9 
and the "thin-layer" solution, Eq. (|22|) , is an approximate 
solution for larger 9. For power-law like spatial correlation 



functions, ^ cx 



a sensible definition for the position 



of the break is at the angle 6'brcak (in RAD) where both 
approximations intersect: 



A, 9 



1-7 



break 



e(%/2v/l -COs6lbrc 



(C.l) 



Using (fTS]) and pT]) one obtains after some algebra: 



'1 - cos( 



break 



r(7/2) 



0Fr(7/2-l/2) 



dr p^{r)r 
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(C.2) 



For 7 S [1,2] the left-hand side of the previous equation 

^brcak • (C3) 



is for ^broak ^ 60° within 10% accuracy just 



nl-7 
break 



'2Vl — COS^break 



